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Flash-algorithm track-reconstruction routines with speed 
factors 3000-4000 in excess those of traditional iterative rou- 
tines are presented. The methods were successfully tested in 
the alignment of the Test Beam setup for the ATLAS Pixel 
Detector MCM-D modules, yielding a 60 fold increase in align- 
ment resolution over iterative routines, for the same amount of 
alocated CPU time. 



I. INTRODUCTION 

In many particle physics experiments high-precision 
alignment can be performed provided there is sufficient 
data and enough computer (CPU) time alocated. While 
statistics may suffice for many experiments, CPU time 
can only be decreased through the use of flash recon- 
struction routines. For instance in pixel detectors, align- 
ment parameters can be tuned iteratively to satisfy good 
track reconstruction in all events. However, reconstruc- 
tion of all tracks in all events, for each iteration in align- 
ment parameter space|^ is CPU exhausting and in prac- 
tice can only be applied to blocks of tracks. If the track 
fits themselves are also iterative double- nesting of 
iteration loops occurs and CPU times reach on the or- 
der of months. On the other hand if track fits are 
semi/analytical, CPU times drop 3-4 orders in magni- 
tude and such an approach becomes feasable. 

The paper starts with the basic fit, examines the va- 
lidity of the standard quadratic approximation for and 
presents analytical and semi-analytical track reconstruc- 
tion algorithms, together with CPU-clocked examples^. 
The methods were developed gradually from the simple 
case of straight tracks (magnetic field free environment, 
or high momentum tracks), to the more demanding case 
of helical tracks. The Test Beam setup of the ATLAS 
Pixel Detector Q MCM-D modules Q was used as pro- 
totype example for the applied methods, although the 
latter methods are general and applicable to a number of 
detectors using tracks. 



II. GENERAL FITS 

The general expression for in a fit is the sum of the 
squared normalised residuals over the set of experimental 
points: 
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where the error of the j -residual ctj;^, is related to the 
direction of the residual itself, Aj through the covariance 



matrix cr" 
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More often however, a quadratic]^ approximation of 



xLact is used: 
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with the equivalent residual error: 



J(A) 



4 • ^7' Ai 



(5) 



If a track impacts a point's error ellipsoid at rt = 
tgifiincid) with respect to one of the principal axes, then 
approximating |Aj) as perpendicular to the trajectory, 
the two cr's can be written as: 



^An example of an analytical alignment procedure {vs. an 
iterative one) is illustrated in |l| for the SLD End-Cap 
Cerenkov Ring Imaging Detector. 

^The tests were performed on a DEC- ALPHA 878 machine, 
with EV 5.6 processor at 433 MHz, 640 MB RAM, and run- 
ning under OSFl V4.0. Coding was performed in FORTRAN. 
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A J ■ A J = Al/al + Al/cTy + Al/cr;, for a] assumed 
diagonal. 
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both reducing to cr^^ or a1 for tracks impacting afong 
one of the principal axes. 

The Test Beam stand for the ATLAS Pixel MCM-D 
modules is a Telescope setup, with 4 tracking elements 
(Sirocco strip Detectors) and two slots for Pixel Module 
evaluation. The setup was mounted on a marble optical- 
bench, itself placed on a rail that allowed it to be moved 
into the active area of a 1.4 Tesla Spectrometer Mag- 
net. The latter was used to determine the magnetic field 
influence (x-shifts) in an environment comparable to the 
ATLAS Detector. The Sirocco strips were mounted along 
Az — 1.4 m of beam- line with a tolerance of ±0.5 mm. 
The strips, 30/im wide, provided a resolution of 4/im in 
the X- and y-directions, while the Pixels, 50/xm x 400/^m, 
a resolution on the order of 14/im x 180/im depending 
on the cast technology of the chips. The error ma- 
trices of the Test Beam stand Telescope points had thus 
associated ellipsoids with aspect ratios of 1:125, the dif- 
ference between fitting a 3D-line to these c's and one to 
spherical cr's being 0.1 //m in the Telescope's mid-plane. 

The difference between Xexact ^^.^ xipprox depends 
strongly on the impact angle of the track onto the in- 
dividual error ellipsoids: 
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For the Telescope setup the track's impact angle was on 
the average 0.15 nirad, with a corresponding Ax^/Xeicact 
on the order of 0.02%. When the tracks impact however 
at an arbitrary angle, Ax^/xLact can reach as high as 
3000 for the current (Jz/cTxy ratio, even if Aa'^ j a^^^^^ is 
on the order of unity. 

For tracks impacting all points at a constant angle 
("stiff' -tracks), Ax^ is constant along the trajectory and 
the two methods yield identical results. If the track is 
measured however piecewise in two different sub-systems, 
or it is composed of an ensemble of points with different 
cr's (different types of detectors), then even for straight 
tracks the two solutions differ. For tracks bending in 
magnetic field, the track's impact angle changes contin- 
uously along the track, Ax^ following as: 
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Within the 1.4 m of the Telescope, the 180 GeV/c tracks 
used bend in the B = 1.4 T magnetic field equivalently to 
Ax' / xlxact — 0.02% up front and 16% downstream, the 
approximative method pulling the fit increasingly tighter 
towards the end - on the order of 0.5 fim per point. The 
effect is evidently insignificant, both in the Telescope 
setup, as well as in the real B-physics context of AT- 
LAS, meaning 2 Tesla magnetic field, tracks of momen- 
tum greater than 1 GeV/c, and a measured track length 
of approximately 0.14 m. Over this span the expected 
point to point change in x' is less than 4 % per %-Ax^. 



III. LINE FITS 

In most cases it is possible to interchange the non- 
linear expression (||) with its quadratic approximation 
(^, allowing analytical solutions to be given for "stiff- 
tracks" i.e. - particle out of magnetic field, weak field 
with respect to track momentum, or distance travelled 
small with respect to existing resolution. 



The simplest fit is for a 



a ■ 1 = const. : 
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Parametrising the tracks as: 



r = fo + Xn (10) 
with fp ~ 1 and fo ■ ft — 0, equation (|^) becomes: 
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where 5j — fj~f\) — Xjn. The minimum condition implies 



locally Xj 



and globally: 

fo = (1 - nn) {r] 



(12) 



where (r) denotes average over measured points, and 
M = (ff) — (r) (r ) the spread ellipsoid of points around 
(r). The 3 eigen- values of M represent the length of 
the track (/io); and the two transversal variances to the 
line-fit. 



For a' — diag{a',(jy,(j1) — const, equation (11) holds 



again, however in normalised form, with — s- cr^^r;, 
To u^^ro and n — s- a^^n. Adapted to the geometry 
of the Telescope this solution has been clocked to 0.033 
/is/2D-fit and 2.6 /is/3D-fit. These times are half or less 
for current 1000-1500 MHz clock machines. CPU-wise 
this allows the alignment to be performed using blocks of 
tracks. The intuitive picture of this procedure is aligning 
two incomplete 3D images: where there is a track, there 
is a pixel on the image. The more tracks in an align- 
ment loop, the more pixels per image, and the better the 
alignment. This idealistic picture is cut cold however, by 
either lack of data, or of CPU power. In most of the con- 
temporary HEP experiments the limiting factor is CPU 
power, as the reconstruction of the objects used in the 
alignment (in this case tracks) can be quite CPU costly. 
Eliminating this problem with fiash-algorithms opens the 
way to high accuracy alignment. The increase in align- 
ment resolution is proportional to \/N uock ■, and the CPU 
cost to Nuock- Holding CPU time fixed, any speed in- 
crease can be equivalenced to a resolution increase. Thus 
a factor of 4000 in speed, would be equivalent to a fac- 
tor of 63 in resolution. To exemplify this, consider the 
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Test Beam setup, where the tracks impact the Sirocco 
planes at almost normal incidence. At precisely normal 
incidence any z-misalignment would be un-noticed. The 
tracks do impact however at an angle on the order of 
0.15 mrad, giving an alignment "lever arm" for Az mis- 
alignments on the order of 0.15 /im/mm. With the exist- 
ing Sirocco resolution (~ 4 fim) the resolving power (per 
track) for Az would be on the order of 27 mm. Increasing 
the number of tracks per one Az misalignment loop the 
resolution improves dramatically, as illustrated in figure 
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FIG. 1. Translational alignment parameter (Az) his- 
tograms for blocks of 2, 4, 8, 20, 40 and 80 tracks, for Sirocco 
Plane 1. The effective number of track-fits in each histogram 
is ~ 3 million. Although in the Telescope setup the alignment 
lever arm for Az is very low (0.15 /xm/mm, with a a = 4 /im 
detector resolution), the plots show a dramatic improvement 
in Az resolution with increasing number of tracks [Nuock) per 
iterative alignment loop. The electronics noise is cut-down by 
a factor of 1/ uock ■ The approach demands however exten- 
sive CPU power, analytical solutions being needed in order to 
keep the problem within the capabilities of existing resources. 
All figures are in mm. 

Any multi-tile pixel detector benefiting from this high 
precision alignment method would better perform in 
physics involving the resolution of vertices. With the 



increasing energy frontier (LHC, TESLA/NLC) particle 
ID techniques are replaced with other means of signal 
identification. Signal "concentrators" that can serve this 
purpose are the B-mesons. Their identification depends 
crucially on vertexing resolution|^ 

An improvement to the example considered would be 
to include in the track fit also the Pixel Demonstrator 
points. This would mean fitting to points with different 
error ellipsoids and the imposibility of "absorbing" all 
ct's into ro and n in an unique way, as done previously. 
For such al ^ (t| cases, the solution is given by a set of 
self-consistent equations: 
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solvable semi-analytically in approximately 3 iterations, 
starting from the previous analytically exact solution. 

The approach used in this paper was to tune the align- 
ment^ by looping over the reconstruction of 2000 tracks, 
the alignment residuals of the Telescope Sirocco plane- 1 
being shown in figure^ (top). The Pixel Detectors under 
evaluation were to first order aligned analytically (exact 
solution), and subsequently tuned in a fashion similar 
to that of the Sirocco planes, in order to compensate 
for the effect of the low y-resolution on the better x- 
resolution. The difference between the Pixel demonstra- 
tor measured and Telescope predicted positions of the 
track hits is shown in figure ^ (middle and bottom) for 
a set of two Pixel Detectors under test. The identical 
szrfe-resolution of the pixels in the x and y-directions (fit 



*For example selecting B-events with high purities and effi- 
ciency can be achieved by applying a minimally missing p± 



correction to the m,r evaluated vertex mass In the AT- 

LAS context pj_ would most likely be referenced to the axis 
of the jet containing the B sub-jet. The method depends 
crucially on the vertexing accuracy, yielding for instance at 
SLD 1^ B-events with sample purities on the order of 91-99 
% and 65-20% corresponding efficiencies. 

^The alignment consisted of two separate stages, firstly out 
of magnetic field, and secondly in magnetic field. In both 
cases the translational and rotational degrees of freedom were 
firstly determined for the 4 Sirocco planes. The Pixel Modules 
were aligned with respect to the Telescope's Sirocco planes. 
The alignment in magnetic field required an external marker, 
which was the setup's trigger: a PIN Diode. Its change in 
image with magnetic field, on a "focal" plane, determined 
the absolute shifts in the tracking elements of the Telescope 
setup (the 4 Sirocco planes), and by reference, those of the 
Pixel Modules under evaluation. 
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parameter P3) - expected from uniform technology on 
the chip - gives credit that the ahgnment conducted to 
physically tangible results. 
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FIG. 2. Final alignment residuals for the Telescope's 
Sirocco Plane-1 (top), Ua; ~ 3 pm, cry ~ 4 /xm, and for 
two Pixel Detectors under evaluation, cast in different tech- 
nologies (middle and bottom). The equal side-resolution 
of the Pixels (fit parameter P3, aside — llfim, essen- 
tially a measure of charge division resolution between pix- 
els) in the x and y-directions gives confidence that the align- 
ment conducts to physically meaningful results. Fit param- 
eter P4 gives the single-hit resolution of the pixel "hat" 
{ahat = 1 ... 10/im, basically the pixel's geometrical resolu- 
tion minus the side-resolution) . The center "thinning" of the 
Sirocco a;-residuals (from a gaussian) by ~ 1 was due to 
the residual magnetic field of the Spectrometer Magnet, as 
Lorentz drift of the charge carriers in the strips. All figures 
are in /im. 



IV. HELIX FITS 

Track fits over small arcs of helices are very sensitive 
to the fluctuations of the experimental points. The ra- 
dius of curvature (giving the momentum, or conversely 
the magnetic field) has consequently a large error, due 
to the points moving within their error bars. To reduce 
the errors of such "weak" parameters, fits to collections 



of tracks - experiencing the same conditions (magnetic 
field in this case) are used. Iterative fit routines consum 
however between 90000-300000 fis/3B helix-fit, the pre- 
viously illustrated method (figure [^) being inapplicable 
in this case. On the other hand, although desired, flash 
semi-analytical algorithms for helix-flts are non-trivial to 
derive. Most of contemporary High Energy Physics ex- 
periments however involve large energies and the helical 
tracks span typically less than 1° of an arc (0.2° in the 
case of the Telescope) . It is possible then in such cases to 
perform 3D-helix fits semi-analytically, by perturbatively 
curving a line-fit into a helix-fit. 
The helix equations are: 



cc 

dtp= -^PX B 
dtf— p/ E 



(14) 



where E is the particle's energy, p its momentum and B 
the magnetic field. The solution to equations (|lj) is: 

f=fo + An+^/(A)Fn+^5(A)Gn (15) 



where A = v^t — uRt is a "linear" distance travelled by 
the particle, Co = \e\(?B /E the helical rotation pulsation, 
fi — pq/po the direction of engagement of the particle 
onto the magnetic field, ub — eB/\e\B, R — pQ/\e\B a 
parameter related to the radius of curvature of the helix 
Rheiix = R\Jl - [n ■ ub)'^, /(C) and g(C) two functions: 
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ff(C) = |(C-sznC) 1 
respectively F and G two tensors: 



(16) 



F = X C 
G^CC 



■ 1 



(17) 



that satisfy Ft = -F, = G, F G = G F = -C^F, 
F2 = G, and G2 = -C^G. The vector C is ub/R. 

It is evident that for i? — s- oo (or equivalently A — > 0), 
expression (|l5| ) reduces to the parametrisation of the line 
( |lO| ) used in performing line fits - which is the require- 
ment for the perturbative approach. In most experiments 
the third order approximation /(C) — 1 and g{(^) ~ 1 
holds up to the following limiting factors: 

• geometric - the arc of helix should not exceed a 
length beyond the approximation validity for /(C) 
and (/(C)- This is related to the demanded resolu- 
tion (T and the particle's momentum: 
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where in the above, A and a are expressed in [m] 
and p in [GeV/c]. The value for the momentum is 
for the Telescope setup. 

dE/dx - the loss of energy along the trajectory 
determines a "tighter" helix, the deviation: 



y-E /dE\ 



(19) 



needing to be smaller than 4 times the allowed tol- 
erance in the Pixel plane. 

multiple scattering - multiple deviations from 
the direction of flight add up to a displacement of: 



a ~ 0.6A6'^ 



(20) 



where A is expressed in [mm] , 9rms in [mrad] and 
a in [/im]. This should determine an error in the 
Pixel plane no larger than the allowed tolerance. 

The semi-analytical helix fit procedure has 3 steps: 

1. Estimation of n, the engagement direction of the 
particle onto the magnetic field. This is obtained 
with small CPU demand via a 3D line flash-fit to 
the first 3-4 points of the trajectory. The vector 
n is an eigen- vector of M = (rr) — (r)(r), hence 
any perturbation (5M = M heiix — M ^ne changes it 
only to second order, and for numerical purposes n 
can be considered constant. 

2. Using the ft found above, the second order term 
corrections to rg and A^ can be estimated: 

AA., = ^{fi- fo) -nxuB 
K 

Afo = (f ) - fo - (A)n - x ub (21) 

ZK 

computations again only modestly CPU demand- 
ing. 

3. Introducing the third order term and using the pre- 
viously corrected (n, ro, Ai), local and global equa- 
tions for the parameters can be written: 

((A-l-yF-f^G)p)=0 

(p)=0 (22) 



Expanding to first order, the corresponding correc- 
tions (An, Afo, AAi) must satisfy: 

a^AX^ + a^■ Afi + h- Afo+ l3i =0 
{a AX) + (A^)! • Af + D Afo + = 
(bAX) +D''An+l- Afo + a ^0 (24) 

where: 

cki = f ^ — f • F f + fo • F n 

Pi — ■ n — fi ■ ft + Xifi^ + Aifo • F ft— 

1 9 1 9 

Xin ■Fn+ -Xfro -Gn- -Xff, ■ G n+ 
^Aff -Gn 

= fo - f + 2XS + XiF fi - XiF fo 
bi = ft + XiF ft 

6 = (A)fo - (Af ) + (A2)f + (AV>- 
l{X^)Fro + \{>?)Gro-\G{X^r) + 

12^ ' 

a = fo - (f) + (A)n + i(A2)Ff + \{X^)Gn 
2 D 

1 



D = (A)1-^(A2)F 



(25) 



By eliminating AAi = — (A + fliAf -I- &iAfo)/ai 
equations (|2j) become: 



M Af + N Afo = f 
N^Af + R Afo = f 



(26) 



where: 



M = (^)-(A2)1 



N==(-)-(A>l + l(A^>F 

,bb. 
R= - -1 
a 



(27) 



and : 



where pi are the residuals of the points to the fitted 
curve: 

A^ A^ 

Pi = -fi + fQ + Xifi+^Ffi+-^Gfi (23) 



a 

TT = (T- ( } 

^ a ' 



(28) 
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To zero* order the three M, N and R tensors are 
proportional to (1— nn), the non-invertable perpen- 
dicular projection to n, by factors of (A^), (A) and 
1. Therefore in the numerical approach, the inver- 
sion is obtained by decomposing the operators into 
a part proportional to (1 — nn) and a "remainder" : 

R = (1 - nil) ■ {2trR - Rnn - Rts)/4 + - (29) 



data reconstruction of the Test Beam results for the AT- 
LAS Pixel Detector MCM-D modules, and can be of 
great impact for high precision alignment and reconstruc- 
tion in all pixel detectors. Such routines could be used 
for instance by the ATLAS Inner Detector to specify 
precision p± corrections to the to^ evaluated mass of 
B-vertices, in order to strongly suppress c-decay back- 
grounds. 



The solution (An, Aro, AA;) is therefore: 
An = (RN-^M - Nt)-i(RN-V- tt) 
At^o = N-i(f- MAn) 



AA,; = {Pi + Si ■ n + bi ■ Afo) 

a,. 



(30) 



All quantities in this section were considered nor- 
malised - i.e. a^^r — ^ r, although for notation 
simplicity they were written as the quantities them- 
selves. 

The CPU demand of the above 3 steps is under 15 
fj,s. For better precision however, the last step can be 
repeated twice, bringing the 3D-hclix fit to 22 /is. This 
is at least 4000 times faster than any iterative version of 
the fit. 

To complete the fit in all generality, dE/dx can be in- 
corporated by perturbatively bending a pure helical track 
into a o?i?/c?a;-helix. This would allow in principle to in- 
fer the track's mass, although with large errors in certain 
energy ranges. 

Using the fit over blocks of 2000 tracks the alignment 
in magnetic field was checked and adjusted. The fit 
was then used on blocks of 10 tracks for a fine scan 
of the beam energy delivered by the SPS to the Test 
Beam setup. This was expressed in normal-impact ra- 
dius of curvature equivalent^, and it is shown in figure ^ 
(bottom-right). 



V. CONCLUSIONS 

Analytical methods were shown to have a dramatic 
impact on the speed of line and helix fits, bringing down 
CPU usage by 3-4 orders of magnitude. The methods 
developed were successfully tested in the alignment and 
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FIG. 3. Normal-impact equivalent radius of curvature for 
single tracks (top-left) and corresponding of fit vs. radius 
of curvature (bottom-left). The resolution was on the order 
of 16 m and the average radius of curvature 310 m. Using 
blocks of 10 tracks (top- right) the resolution improves beyond 
the statistical factor, to ~ 3.3 m, the stability of the SPS 
delivered beam being tracked vs. time bottom-right. 
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